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Abstract 

We suggest and discuss a simple model of an ideal gas under the piston to gain an insight into the workings of 
the Jarzynski identity connecting the average exponential of the work over the non-equilibrium trajectories with 
the equilibrium free energy. We show that the Jarzynski identity is valid for our system due to the very rapid 
molecules belonging to the tail of the Maxwell distribution. For the most interesting extreme, when the system 
volume is large, while the piston is moving with large speed (compared to thermal velocity) for a very short time, 
the necessary number of independent experimental runs to obtain a reasonable approximation for the free energy 
from averaging the non-equilibrium work grows exponentially with the system size. 

1 Introduction 

The celebrated Jarzynski identity is perhaps the most recently discovered simple general formula in elementary 
statistical mechanics: 

(e w/kBT ) = e ~ AF/k * T (1) 

(see below about sign convention). The claim of this charming simple formula is as follows. Suppose we have 
an arbitrary system and let us consider two states of this system specified by parameters, say j4i n i t i a i and ^4a n ai; 
these could be volumes, magnetic fields, or just about anything else. If the system comes to thermodynamic 
equilibrium at A inMa i, it has free energy -F(Ai nitia i) = F init i a i; if the system is equilibrated at Afi na i, its free energy 
is -F(j4fl na i) = -Fflnai. The difference of these free energies is AF = Jfi na i — Finitiai- According to elementary 
thermodynamics, if we drive the system from initial to final state by a reversible process, such that the system 
remains at equilibrium at every stage, then we have to perform work, —W, which is equal the free energy change: 
—W = AF; if, on the other hand, the process is not reversible, then the second law of thermodynamics tells us that 
(-W) > AF. We know, of course, that the second law of thermodynamics is of a statistical nature, and, therefore, 
from time to time, very infrequently, the fluctuations occur in which — W < AF. These fluctuations might be very 
rare, but with large W (strongly negative —W) their contribution to the average of e w l k B T might be significant. 
The Jarzynski formula Q tells us precisely that: when all fluctuations, including those violating — W > AF, are 
taken into account, then the average of (e w ^ kBT ) reduces to e _AF//fe - BT . 

It is worth repeating that — W is the work performed by an external force on the system; in other words, W is 
the work performed by the system itself. We use this sign convention (perhaps somewhat non-standard) because it 
will make for positive W and save us some writing in interesting cases below. 

The tempting use of this result is to circumvent in computer simulation or in real experiment the often painful 
stage of equilibrating the system. Instead, it should be possible to run the system many times without any worry 
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of the equilibrium, repeatedly measure the work, W, and still obtain the equilibrium information, AF from the 
formula Q. The problem is that it requires exploration of all sorts of fluctuations along the way and measurement 
of the work W for every fluctuation. 

There are already quite a number of works exploring various aspects of the Jarzynski formula. Papers [I] 
address mathematical foundations of the Jarzynski identity in the context of fundamental statistical mechanics. 
Some authors describe the Jarzynski identity in terms of 'transient violations of the second law of thermodynamics' 
6 . There is a history of cases when such fluctuation effects were mentioned under the name of 'temporary violations 
of the second law' [7||Sj, but many people feel that it is incorrect terminologically to say that fluctuations violate the 
second law, even if temporarily. Here, we do not make any firm commitment to any terminology in this sense, and 
just say that the Jarzynski identity is based on proper exploration of a representative set of fluctuations. Various 
ways to apply the Jarzynski relation in computer simulations and in experiments are discussed in |o1 1§1 I1L)( 1111 
(see also extensive literature cited in these works). 

Nevertheless, we found that one aspect is missing in the current literature, namely, the pedagogical aspect. The 
Jarzynski formula appears so simple and so general that there ought to be a simple way to explain it and to gain 
an intuitive insight of it on a very elementary level. This motivated us to look for the simplest possible example in 
which the Jarzynski equation shows some non-trivial results. Since the conceptually simplest subject in statistical 
mechanics is undoubtedly the classical ideal gas, the goal of the present paper is to work out the application of the 
Jarzynski formula to the ideal gas. Of course, we shall find nothing really new in terms of factual results, but we 
hope for a good new insight. 

2 Model and formulation of the seeming paradox 

Consider some amount of ideal gas in a vessel under a piston. Everything is supposed to be in accord with elementary 
physics textbooks: ideal thermoisolation, mass-less piston moving without any friction, etc. Suppose initially the 
piston is some distance L from the bottom of the vessel, and that the gas temperature is T. Let us now move the 
piston by some distance AL and stop it again, thus preparing the final state. 

The following way of thinking seems quite logical from a physicist point of view. To make the situation dramatic, 
let us suppose that we move the piston at a very high speed, much faster than the speed of sound in the gas, or, 
in other words, much faster than the averaged thermal velocity of the molecules. Then we can roughly say that no 
molecules will be able to chase the piston while it is moving, no molecules will hit it, and there will be, therefore, 
no work. When the piston is stopped at the end, molecules start arriving, they do bombard the piston, but since 
the piston does not move at this stage, the work is still zero. This logic leads to the conclusion that in this case 
W = 0, implying (e w/kBT ) = 1, while obviously AF / 0, which seems to contradict the Jarzynski identity Q. 

To resolve this paradox we have to remember about the tail of the Maxwell distribution: however large is the 
speed of a piston, there is still some probability of molecules moving fast enough to chase the piston and hit it 
while it is moving. This already suggests that the Jarzynski identity Q has to do with the tails of the relevant 
distributions. To make this statement more precise, we shall compute the probability distribution of the work W 
for our elementary model. This is obviously much more than just computing the average involved in Jarzynski 
formula Q. 

In order to make our article more pedagogical, we shall start with proving the very identity (0 for our specific 
system. We shall also relegate cumbersome calculations to the appendix. 
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3 Calculations 



3.1 Average value of e^ fcflT 

Since we plan to consider an ideal gas, all molecules will contribute to both W and AF independently. Therefore, 
we can imagine the Jarzynski formula Q re-written as 

((e w ))" = ( e - AFl/kBT ) N , (2) 

where W\ and AF\ are the work and the free energy change per one molecule. We see that for the ideal gas, 
quantities W\ and AFi satisfy the Jarzynski formula looking identical to fl). We, therefore, restrict ourselves for 
simplicity to the "ideal gas" of just one molecule, and also for simplicity we suppress the index 1 in writing W and 
AF. Thus, we keep considering formula |{T). but we think now about just one molecule in an ideal gas. 

Furthermore, to simplify writing, we assume that the temperature is such that ksT = 1, the mass of the 
molecule is m = 1, and the piston is moving during the time interval r = 1. 

Figure illustrates the system consisting of a thermally-isolated cylinder, a piston moving at speed v p , and a 
single molecule initially at position x with velocity v. The molecule bounces off the walls elastically, so we are 
concerned only with the one-dimensional motion indicated. The space-time diagram depicts the trajectory of the 
piston and the trajectories for a molecule initially moving toward the piston (dashed lines) and for a molecule 
initially moving away from the piston (thinner dashed lines). 

Let us focus on the work W done by the single molecule on the piston, in a time interval r = 1. The quantity e w 
(fesT = 1) is to be averaged over the possible initial states (a;, v) of the molecule drawn from a Maxwell-Boltzmann 
distribution. Therefore, 

(e w ) = ^ =^2£ (3) 

ff dx f°° dve-v 2 /? 

JO J — oo 

where w T (x,v) is the work done by the gas molecule given the initial coordinate (x,v) and time elapsed r. To 
compute this function w T (x,v), we first need to work out the collision times between the molecule and the piston 
and the work done after n collisions. 

Let us first assume a positive initial velocity, in which the molecule can strike the piston first before hitting the 
left end of the cylinder. The time taken for the first collision with the piston is ii = _ g . After the collision, the 
velocity of the molecule relative to the piston gets reversed and the speed of the molecule gets diminished to v — 2v p 
(assuming v > 2v p ). The time taken for the second collision with the piston is given by = ^rf^f ■ in general, for 
the n th collision 

,+ _ ( 2n -l)L-x 

v- (2n - l)v p ( ' 

Similarly, for a molecule with a negative initial velocity, 

_ (2n - 1)L + x 
n v- (2ti - l)v p { ' 

These relations can be inverted to give conditions that should be satisfied by the speed of the molecule in order 
to result in exactly n collisions with the piston within a time interval r = 1. For positive initial velocities, 

(2n - 1)(L + v p ) - x < \v\ < (2n + 1)(L +v p )-x (6) 

For negative initial velocities, 

(2n - 1)(L + v p ) + x < \v\ < (2n + 1)(L + v p ) + x (7) 

The work done by the piston on the molecule after one collision is the change in momentum of the molecule 
times the velocity of the piston, 

— Wi = (— (t! - 2v p ) - v)v p = -2(v - Vp)v p (8) 
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In general, the work done after n collisions is 

— w n — —2vv v n + 2Vpti 2 (9) 
Note that the work done can also be calculated from the change in kinetic energy after n collisions, 

- w n = -(v - 2nv p ) 2 - v 2 /2 = -2nvpv + 2n 2 vl (10) 

The work done by the molecule on the piston is positive for an expanding volume. 

We are now facing the laborious task of calculating the integral in the numerator of formula J3J . It is cumbersome, 
because it must include the summation over all possible numbers of bounces of our molecule from the piston. Note 
that a large number of bounces correspond to a very large initial velocity of the molecule, as it has to have time 
to chase the piston for n bounces, even though it looses momentum and gets slower at every bounce. The actual 
calculation is described in appendix^ Using the simplified result 13U1 . the sought average is 



(e W ) 

r dx | due— '!* 

JO J — oo 
T _L ni 

(11) 



r +v * dx f°° dye-" I 2 

JO J — oo 

f L dx f x dve- 2 ^ 

JO J — oo 

L + Vp 



L 

which can be recognized as the ratio of the partition functions at the final (after time r = 1) and initial volumes at 
the initial temperature T = Z(L + v p t,T)/Z(L,T). This expression is also identical to that obtained from 

the Jarzynski identity. 

3.2 Probability distribution of the work W 

The prescription for evaluating the distribution is 

-. t'L /*oo 

P(W) = — / dx I dve- v2/2 S(W -w T (x,v)) . (12) 

V2tyL Jr. J _^ 



The calculations are presented in detail in appendix|B| With the expression Ijlfijl for n, the number of bounces, one 
can get rid of the summation over this number and the distribution function simplifies to, 

P(W) = S(W)P + = —f(W). (13) 



Here, Po is the probability to obtain vanishing work because the molecule is unable to chase the piston or hit it 
even once, 

Po = -jL^ / dx / dve~ {v - x) 12 , (14) 

V2ttL J J_ (L+Vp) 

and the function f{W), which we call the overlap factor, can be formulated as follows (see also figure 0: 

f + when (n-l)(v p +2L)<^-<(n-l)(v p + 2L)+2L 

f(W)=i 1 when (n- l)(v p + 2L) + 2L < < (n- l)(v p + 2L) + 2L + 2v v 

I (i+l) + 1)- 4^T when (« ~ 1)(«p + 2L) + 2L + 2v p < ^ < (n + l)(v p + 2L) 

(15) 

Here, the integer n (which is the number of bounces by the molecule against the piston) is obtained in appendix 151 
and is given by the formula 



I I — 

1 v p (2L + Vp) 



(16) 
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where [. . .] means integer part of . . .. For example, simple algebra indicates that as long as W < 4v p (2L + v p ), we 
have just one collision, n = 1. For the values of work W in the next interval, 4n p (2L + v p ) < W < Y2v p {2L + v p ), 
we have n = 2, etc. 

Thus, the probability distribution P(W) consists of a 5-function peak at W = and a tail at positive W. 



3.3 Limit of large volume and fast moving piston 

As we said in the beginning, the most interesting case is when the piston moves fast, such that hardly any molecule 
can chase it and produce non-zero work. That means, the Jarzynski identity in this case relies exclusively on the 
far tail of the Maxwell distribution. Let us consider the probability distribution P(W) in this limit, v p ^> 1. 

It is reasonable to assume simultaneously that the volume is large enough, such that L 3> v p . Since in more 
traditional units this condition reads L 2> « p t, it means that the piston moves fast, but for a very short time. 

In this case, the distribution is dominated by the single bounce, that is, n = 1. Assuming n = 1 and v p <C L, 
we have 

/ when 0<W < Av p L 

f(W) ~ | w_ when 4VpL<W<8VpL ( 17 ) 

which yields the probability distribution expression 

e -k( v p+^) w 

P(W) ~ S(W)P + = -£-. (18) 



valid at i > j)p > 1 and W < 4Lw p . For the larger W tail of the distribution P(W), we have to include the second 
line of eq. 1171 . and for even larger W also higher n values. Luckily, there is no need to do that, because the simple 
approximation I18H is good enough to capture the Jarzynski result. Although the Jarzynski identity involves (e w ), 
which includes integration over all values of W, the integral converges rapidly enough to yield the correct answer 
within the region of applicability of formula I18H : 

P(W)e w dW = P + -^L- e ^ " x Jl-dW 



L dx 1 
o L 



+ 



1 r°° _i(w-2v 2 ) 2 2v 2 

e * x - — --dW 



= 1+^. 

Here, we made approximations in both the Po term, by extending the integral limits to (— oo,oo), and in the tail 
term, by setting the integral lower limit to — oo. This way, we do recover the Jarzynski formula J I I i . 

Let us also calculate the probability of obtaining non-zero work values (one or more collisions) as well as the 
average work done, i.e. the th (without the Po term) and I s * moments of the distribution. Using expression 1181 . 
we have 



o+ 



l\V KrW 2 

P{W)dW = / / _ r - / We 8v p dW (19) 
(W) = / WP(W)dW = / / _ r - I W 2 e 8v p dW (20) 



n w2tyLv p j q 
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For large v p , we neglect the term W /(8v p ) in the exponents upon integration, yielding 



P w>0 ~ —}—e~^ (21) 

V 2tvLv p 

(W) ~ ^—e-^ 2 (4)=4P w>0 (22) 

V 2ivLvp 

To the same approximation, the probability to obtain zero work W is equal to 
As one could have expected, these probabilities are governed by the tail of the Maxwell distribution. 

3.4 Comparison with simulations 

Measurements were made in computer simulations and compared with the results obtained in the previous sections. 
The conditions for the 'pulling' experiments were as follows: the 'pulling' time r = 1, the temperature was set to 
ksT — 1 (which sets the width of the gaussian distribution from which the initial velocities were selected), and the 
number of trials or iterations used per measured average was 100, 000. The parameters that we varied were the 
piston velocity v p and the initial piston length (or 'volume') L. 

Figures [3] and [I] present the distribution of probabilities P(W)dW for two different sets of piston velocities and 
piston lengths. In the first case, the piston velocity was set to v p = 0.01 and in the second case was set to v p = 1. 
The effect of the overlap factor is evident in the former case of a slower moving piston. 

Figure presents data for the average (exp(WyfcsT)}. The Jarzynski identity predicts that this average should 
be (1 + Vp/L) = 2. 

Figure [(J presents data for the average work done. The expression corresponding to a free energy change at 
constant temperature T, namely In (1 + v p /L) = ln(e ) = — AF, is plotted, as well as the expression for the average 
work done 1221 for large v p (v p 3> 1) . As the velocity v p increases, the average work done is seen to shift from 
one regime (in which (W) ~ — AF) to another (in which (W) < —AF). If we take — (W) — AF = Wa; s as some 
measure of 'dissipation', then it is also seen that this quantity increases as v p increases, although the difference is 
not much more than ksT . 

In figure [7| the expression for the probability of obtaining non-zero work values in the high velocity limit is 
compared with the fraction of trials in which a collision occurred between molecule and piston. Due to the rarity 
of collisions (and dominance of single collisions) in the range of velocities tested (v p — 1, 1.5, 2, 2.5, 3, L — 1), the 
fraction of trials with collisions is identical to the average number of collisions. 



4 Discussion and conclusion 

Let us look closer at our main results obtained in section l!01 One question to ask is this: how many times should 
one perform the experiment of moving the piston in order to get a reasonable estimate of the average {e w )7 At the 
very least, in order to get the non-zero answer for the free energy difference from Jarzynski formula Q, one has to 
get at least one case of non-zero work. For this, one has to perform about l/P\y>o experiments, which is already a 
very large number at v p 3> 1. In fact, as our calculations show, in order to recover the Jarzynski identity, we have 
to continue integration into the region where W is as large as about Lv p . In practical terms, this means, we have 
to perform as many runs on the system as to get at least a few realizations with the work of this order. According 
to the formula 118fl . the corresponding probability is roughly proportional to e~ 2L . In other words, this requires 
about e +2L runs. Restoring the more traditional notations with fcsT and r, we estimate the necessary number of 
runs (or trials) as exp [mL 2 /r 2 k B T]. Clearly, this is a very large number. 
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In practical terms, one may also want to know if the use of the Jarzynski identity is useful. At first glance, 
it seems extremely useful: one apparently does not have to equilibrate the system and by doing purely non- 
equilibrium measurements, one nevertheless recovers the equilibrium free energy. Our example suggests that the 
situation might be a little more tricky. Indeed, to do equilibrium measurements, one has to proceed very slowly, 
to keep the system close to equilibrium all the time; for this, r has to be larger than the system relaxation time, 
which grows with the system size L (in our dimensionless variables, this corresponds to the limit » p < I C 1). 
But, on the other hand, if one proceeds very rapidly, then one has to perform exponentially many experiments in 
order to catch the exponentially rare but decisively important fluctuations. This consideration suggests that there 
might be some optimal strategy. For the ideal gas model, such optimal strategy is most likely the (classical) slow 
'equilibrium' experiment, because the time for such an experiment grows only linearly with L, while the time for 
a 'fast' experiment is exponential. For other systems, the optimal strategy might be intermediate between one 
very slow experiment and very many rapid ones. Unfortunately, it is clear that such optimal strategy is highly 
sensitive to the particularities of the system in question. Our model, which is an ideal gas, is a system with a 
flat energy landscape. For other energy landscapes one may wonder about the trajectories visiting various valleys. 
Unfortunately, the knowledge of these valleys is exactly what one wants to learn from making measurements of 
equilibrium free energies. In any case, the insight we can gain from our primitive model is that rapid non-equilibrium 
measurements are not automatically advantageous. 

To conclude, we have presented a very naive simple model to look at the Jarzynski identity. We do recover 
the identity for our model, and we are able to demonstrate that its validity relies on the far tail of the Maxwell 
distribution, in the sense that the dominant contribution is provided by the very rapidly moving molecules. We are 
also able to estimate how many independent experimental runs are necessary to obtain the equilibrium free energy 
from the Jarzynski identity with a reasonable accuracy, this necessary number of trials appears exponential in the 
system size. 
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Using the expression for the work done in n collisions 1101 and in view of the inequalities JSJ and Q developed in 
section UTTl we have 



A Computing the numerator in Equation [3] 



i 





(2n+l)(L+v p )-x 



^v 2 e -(-2vv p n+2v^n 2 ) 



dve 



=1 J {2n-\)(L + v p )-x 
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r( L + v p)~ 
l dx I 

JO J-(L+v 

/ ^ U 

Jo n=l [ J C2n-l 

r(2n+l 
J (2n-l) 

j dx j 

JO J -{L + Vp 



dve-i v e° 



-(L+v p )- 

r(2n+l)(L+vp)- 

dve-^ v ~ 2vpn) 

l)(L+v p )—x 

r(2n+l)(L+v p )+x 



dve -l (v -2 V pnf I + 
i-l)(L+v p )+x J 
-L p(L+v p )-x 

dve-* v 

-(L + v p ) — x 



Employing a change of variable, v' = v — 2nv p , 

f (2n-\-\)L — x-\-Vp 



JO „ =1 [J(2 n -1) 

I' 

J (2: 



dv'e-2 V + 

)L — x — v p 
(2n+l)L+x+v p i /2 

dv'e~^ v 

(2n — l)L-\-x — v p 
L pL — x + v-p 



dx I dve~i v2 

J — (L + x) — Vp 

(24) 

If v p is zero or absent, the separate integrals could be coalesced into a single integral of a gaussian from — oo to oo 
and the result is trivial (identical to the denominator in the average). Therefore, let us separate the 'excess' from 
the trivial result, 

/* L con 

, -it, 2 

dve 2 + 



°° ( p(2n+l)L-x+Vp A2n-l)L-x 

^2 I dv+ dv+ 

n=1 I J (2n+l)L-x J (2n-l)L-x 

r(2n+l)L+x+v p p(2n-l)L + x 

/ dv + 

J(2n+l)L + x J(2n 



dv \ + 



' (2n+l)L + x J (2n-l)L + x-v p 



L~x-\-v v pL-\-x-\-v p 

dv+ dv} e~^ v 

L—x J L+x 



The first, third, fifth and sixth integrals after the summation symbol (those with upper limits "... + v p ") can be 
combined after making the change of variables, x' = (2n + 1)L — x, x' = (2n + 1)L + x, x' = L — x, x' — L + x, 
yielding 

r°° rx'+v p i 2 

h = / dx' / dve~^ v (25) 

Jo J x' 

Combining the second and fourth integrals after the summation symbol (those with lower limits "... — v p ") in a 
similar away, after the change of variables x' — (2n — 1)L — x, x' = (2n — 1)L + x, the result is 

h = I dx' I dve'^ 2 (26) 
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Performing yet another change of variables, v' = v — x\ 

poo rvp 

h = dx' dv'e-^ v ' +x "> (27) 



JO 



h = / dx' dv'e-^ v ' +x ' )2 (28) 



o 



(Notice that the series of variable substitutions effectively exchanged the infinite limits associated with v with the 
finite limits associated with x.) After performing the change of variable for I2, v' = —v", and combining the two 
integrals, 

h + h = I dx' J dv'e-^ x ' +v "> 2 + 

dx' / dv"e-^ x '- W,) 
Jo 

dv' I dx'e~^ x '- v ' )2 

J — 'X 

Exchanging the "roles" of x' and v' , i.e. letting v = x' and x — v' , 

h + h= dx dve~2 v (29) 







Using this simplified result, the numerator in the average becomes 

I = I dx dve~^ v2 + dx I dve~^ y2 







L+v p 

dx / dve~? v ~ (30) 



1 ..2 



B Computing the probability distribution P(W), Eq. ( 1121 ) 

The inequalities iJSJ and (|7J developed in section |3~T1 lead to the following partition of the integral 

-L °° p(2n+l)(L + v p )-x 



P(W) = . / dxy] dve~ v /2 x 5(W - (2vv p n~2vln 2 )) + 

v2"kL J n=1 J( 2n -i)(L+v p )-x 

~j=^ I dx ^2 dve ~ v 12 x s ( w ~ ( 2vy p n - 2v W)) + 

V2ttL J J (2n-l)(L + v p ) + x 

1 



L ML+v p )-x 2 

dx I dve~ v /2 8(W-0) 

J —{L+v p ) — x 



Call the first term I\ and the second term I2 (the third term is 'trivial'). Performing a change of variable to 
remove x from the limits, integrating over x and taking advantage of the delta function results in 



h = y -= x — {overlap between 

^-^ y2irnv p 2L 
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nvp + 



W 
2nvv 



, nv p + 



W 

2nvr, 



and similarly for h, 



h = 



and [(2n - 1)(L + v p ), (2n + l)(L + v p )]} 



x — {overlap between 



71=1 



2irnv v 



W 

2nv p 



L, nv p + 



W 



2nv v 



and [(2n - 1)(£ + v p ), (2n + l)(L + v p )]} 



Ii and h can be combined as follows, 

h + h ~- 



n=l 



2irnv v 



x — — {overlap between 



IT 



L, nv p + 



IF 



2nv v 2nv v 
and [(2n - 1)(L + (2n + 1)(L + u p )]} 



E 



2nnv v 



x /( n ,W) 



where the overlap factor / satisfies < / < 1, since the range of the smaller interval is at most 2L. f is also zero 
for negative W, or positive work values —W done by the piston. 

The conditions that must be satisfied by W in order for the overlap associated with integer n to occur are 



2nv p (2(n - 1)L + (n - l)v p ) < W < 2nv p (2(n + 1)L + (n + l)v p ) 



(31) 



Notice that the left boundary of the interval is a function of n(n — 1), while the right interval is a function of 
n(n + 1). Therefore the right boundary can be transformed into the left-boundary by making the replacement 
n — » n — 1. This implies that the intervals 13111 are contiguous and nonoverlapping and that at most one term in 
the summation in P(W) survives. One can solve for the integer n by taking the integer part (or floor function) of 
a solution to a quadratic equation, 

W > 2nv p (2(n+ 1)L + (n + l)v p ) (32) 
which results in formula HIGH in the main text. 
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Figure 1: System and space-time diagram. Solid line - piston trajectory; Thick dashed line - molecule trajectory with 
positive initial velocity; Thin dashed line - molecule with negative initial velocity 
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Figure 2: The structure of the overlap factor f(W), which modulates the exponential in the distribution function. 
This factor becomes a rapidly oscillating function in the limit of small piston velocities. 
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Figure 3: Simulation results (pluses +) together with theoretical calculation (expression l|13fl . crosses x) for the work 
distribution. For each trial run, the cylinder volume was doubled (v p — 0.01, r = 1, L = 0.01). The bin width used 
was Aw = fc^T/lOOO = 1/1000. The average number of collisions between molecule and piston was about 20. 
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Figure 4: Simulation results (pluses +) together with theoretical calculation (crosses x) for the work distribution. For 
each trial run, the cylinder volume was doubled (v p — 1, t = 1, L = 1). The bin width used was Aw = = 1/10. 

The average number of collisions between molecule and piston was low, about 0.08. 
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Figure 5: Evolution of the average of exp(w/kBT) with the number of trials. 
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Figure 6: Plot of measurements of the work done (pluses +) together with the expression In (1 + v p t/L) (i.e. l —AF\ 
with r = 1 and L — 1) and the expression for the average work done for large v p l|22|l . 
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Figure 7: Plot of the expression (|21Jl fline). the fraction of trials with a collision (pluses +) and the average number of 
collisions (crosses x ) . The latter two coincided. 
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